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A new method to study the effects of neutrino masses on a supernova neutrino signal is 
proposed. The method relies exclusively on the analysis of the full statistics of neutrino 
events, it is independent of astrophysical assumptions, and does not require the observa- 
tion of any additional phenomenon to trace possible delays in the neutrino arrival times. 
A statistics of several thousands of events as could be collected by SuperKamiokandc, 
would allow to explore a neutrino mass range somewhat below 1 eV. 



1. Introduction 

Already long time ago it was realized that Supernova (SN) neutrinos can provide 
valuable informations on the neutrino massesPJ The basic idea relies on the time- 
of-flight delay At that a neutrino of mass m„ and energy E v traveling a distance L 
would suffer with respect to a massless particle: 

At _ 1 / 5.1ms \ / lOMeV \ 2 / m v \ 2 

T ~ v ~ \10k^ J \ E v ) VTeVJ ' ( ) 

where for ultra-relativistic neutrinos 1/v — E v /p v ~ 1 + m^/2i?^ was used. The 
dispersion in the arrival time of about twenty v e from supernova SN1987A was used 
in the past to set the model independent limit m Pc < 30evP while a recent de- 
tailed reanalysis obtained, within the SN delayed explosion scenario, m Pe < 5.7eV"PI 
Since SN1987A, several efforts have been carried out to improve the sensitivity of 
the method, while awaiting for the next Galactic SN explosion. Often, these ap- 
proaches rely on "timing" events related to the collapse of the star core, that are 
used as ben chm arks for measuring the neutrino delays. The emission of gravita- 
tional wavesp^the v e neutronization burstj^the initial steep raise of the neutrino 

luminositypl and the abrupt interruption of the neutrino signal due to a further 
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collapse into a black holt have been used to this aim. However, there are some 
drawbacks to these methods: firstly only neutrinos with arrival time close to the 



*Work done in collaboration with J.I. Zuluaga, and partially supported by Colcicncias under 
contract 1115-05-13809. 



5, 2008 6:33 Proceedings Trim Size: 9.75in x 6.5in 



Main 



2 

benchmarks are used; secondly the observation of the benchmark events is not al- 
ways certain, and in any case some model dependence on the details of the SN 
explosion is generally introduced. A new method that is free from these drawbacks 
was proposed in Ref. El The method relies only on the measurement of the neutri- 
nos energies and arrival times and uses the full statistics of the detected signal. It 
is also remarkably independent of particular astrophysical assumptions since no use 
is made of benchmarks events. The basic idea is the following: in the idealized case 
of vanishing experimental errors in the determination of the neutrino energies and 
arrival times, and assuming an arbitrarily large statistics and a perfectly black-body 
neutrino spectrum, one could use the events with energy above some suitable value 
E* (to suppress the mass effects) to reconstruct very precisely the evolution in time 
of the neutrino flux and spectrum. Once the time dependence of the full signal 
is pinned down, the only parameter left to reconcile the time distribution of the 
low energy neutrinos with the high energy part of the signal would be the neutrino 
mass, that could then be nailed to its true value. Of course, none of the previous 
conditions is actually fulfilled. In Cerenkov detectors as SuperKamiokande (SK) the 
uncertainty in the energy measurement is important and must be properly taken in 
to account. The statistics is large but finite, and finally, the SN neutrino spectrum 
is not perfectly thermal!^ Nevertheless, a good sensitivity to the mass survives, 
allowing to disentangle with a good confidence the case m v — from m v = 1 eV. 

The idea outlined above was tested in Ref. 8 by proceeding in two steps: i) F irst 
a set of synthetic neutrino signals is generated, according to a suitable SN modelP^ 
Neutrinos are then propagated from the SN to the detector assuming two different 
mass values: m u = and 1 eV. ii) The detected signals are then analyzed with the 
aim of disentangling the two cases m„ = and 1 eV. Only the SN-Earth distance 
is assumed to be known, while all other quantities, like the spectral functions and 
the detailed time evolution of the neutrino flux are inferred directly from the data. 

For the time evolution of the neutrino flux and average energy the results of 
the SN explosion simulations given in Ref. E3 were used. The neutrino energy 
distribution was modeled by a Fermi-Dirac (FD) spectrum with time dependent 
spectral temperature T(t) and a 'pinching' factor 77(f) introduced to account for 
spectral distortionJ^. Each neutrino is labeled by its emission time t v and by its 
energy E u , and the corresponding positron produced through the reaction D e p — > 
e + n is also identified by a pair of values (E e , t e ) generated by taking into account 
the detection cross section a{E v ) and the SK energy threshold and resolution. 

The synthetic signals are then analyzed by means of the likelihood function 

C = S(e; T(f), 77(f)) x $(f + St; b, d, /) x a(e) , (2) 

where e is the neutrino energy inferred from the positron energy. The effective 
temperature T(t) and pinching 77(f) of the FD spectral function S, and the detailed 
shape of the parametric function <f>(f + St; b, d, f) that describes the neutrino flux 

Ixl 

are directly derived from the data™ Given a value of m„ the time delay of each 
neutrino is computed according to its energy ej, and subtracted from its arrival 
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Table 1. Results of the fits to the neutrino mass. 
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time ti . For the new array of times — log C is then evaluated, and minimized with 
respect to the flux shape parameters, until the value of the mass corresponding 
to the absolute minimum is found. The power of the method for disentangling 
m„ = 1 eV from m„ = was studied by assuming in turn the two energy thresholds 
E tT = 5 and lOMeV, and two SN-Earth distances L = 10 and 20kpc. For each 
case 40 samples were analyzed. Table 1 summarizes some of the results. The first 
raw gives the percentage of times in which the 95% c.l. lower limit m v is larger than 
the input mass fh v . The second raw refers to the cases when the upper limit m" is 
smaller than m v . Number in parentheses correspond to L = 20kpc. These figures 
characterize the percentage of 'failures' of the method, that therefore appears to be 
reliable in about 90%-95% of the cases. The third row gives the percentage of times 
when m„ = is excluded at 95% c.l. when the signal is generated with rh v = 1 eV. 
We see that for E tI = 5MeV and L = lOkpc the method is successful in more 
than 50% of the cases. The results in table ^ were derived relying on two main 
simplifying approximations: 1) the neutrino energies were generated assuming a 
'pinched' FD spectrum and fitted with a similar two-parameters energy distribution; 
2) no effects of the neutrino oscillations were taken into account. The effects of 
these approximations have been analyzed in Ref. |8] by carrying out 40 simulations 
using the numerical spectra given in Ref. OH and assuming a mixed Ve-v^ composite 
spectrum as would result from neutrino oscillations. It was found that even when 
the u e -v^ spectral differences are assumed to be unrealistically large, the sensitivity 
of the method to the neutrino mass is only mildly affected. 
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